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Ultrasonic control of terahertz radiation via lattice anharmonicity in LiNbOg 
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We propose a novel tunable terahertz (THz) filter using the resonant acousto-optic (RAO) effect. We present 
a design based on a transverse optical (TO) phonon mediated interaction between a coherent acoustic wave 
and the THz field in LiNb03. Wc predict a continuously tunable range of the filter up to 4 THz via the 
variation of the acoustic frequency between 0.1 and 1 GHz. The RAO effect in this case is due to cubic and 
quartic anharmonicities between TO phonons and the acoustic field. The effect of the interference between 
the anharmonicities is also discussed. 



In recent years terahertz (THz) radiation has become a 
tool to study a wide variety of phenomena, necessitating 
the development of experimental techniques to access this 
spectral bandi Radiation in the THz range has been used 
to study picosecond phonon dynamics, polariton propa- 
gation, optical properties of various materials, and has 
been used for imaging purposes There are also sev- 
eral commercial applications that cover a variety of fields, 
including sub-mm wave astronomy, chemical recogni- 
tion and biomedical imaging for disease diagnostics, 
THz imaging and sensing for security applications^— 
As a result, spectrally resolved control of THz radi- 
ation has become an important research topic. One 
implementation uses optical control of carrier densities 
in type-I/type-II GaAs/AlAs multiple quantum wells at 
cryogenic temperatures^ Other implementations include 
magnetically tuned liquid crystals in metallic hole arrays 
and Lyot and Sloe filters, which were shown to be tun- 
able over various ranges between 0.1 and 0.8 THzf^— A 
relative lateral translation of two metallic photonic crys- 
tals has also resulted in a THz filter (tunable between 
0.365 and 0.385 THz) M 

In this Letter we show that the concept of the resonant 
acousto-optic (RAO) effect could be used to produce a 
tunable THz filter. We present a design and model nu- 
merically a LiNbOa-based filter which is tunable in a wide 
range of up to 4 THz. The RAO effect which we use is in 
general a mediation of the interaction between an acous- 
tic wave and a light field by a solid state excitations, when 
the light field and the excitation are in resonance.— ~— In 
ionic crystals with soft TO-phonon modes, the interac- 
tion between a coherent acoustic wave and the THz field 
is mediated by a TO phono m 17 ' 18 The TO phonon is cou- 
pled to the acoustic wave via the anharmonicity present 
in the interatomic potential of the crystal lattice. This 
makes LiNbOa an excellent candidate for the filter, as 
it is both ionic and very strongly anharmonic. In fact, 
the strength of the anharmonicity in LiNbOa is such that 
both the cubic and quartic terms of the Taylor expansion 
of the interatomic potential have to be considered. 

The design of the LiNb03-based tunable filter is shown 
schematically in Fig. 1. A thin slab of LiNbOs (1mm 
thick or less) is placed between doped semiconductor con- 
tacts producing a resonator for ultrasound waves but at 
the same time transparent for THz light. Owing to the 
phonon anharmonicities, two counter-propagating TO- 
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FIG. 1. (Color online) Schematic design of a LiNbOg-based 
THz filter with an acoustic resonator. 



phonon polariton waves in LiNbOs excited by an external 
THz light can strongly couple to each other through the 
coherent acoustic field in the resonator. This happens 
at specific frequencies determined by the resonant Bragg 
condition^ At these resonant frequencies, acoustically 
induced transitions between TO-phonon states result in 
band gaps in the polariton dispersion creating peaks in 
the reflectivity. The anharmonic scattering channels re- 
sponsible for these transitions are "TO phonon ± one 
(two) longitudional acoustic phonon(s)— > TO phonon", 
for the cubic (quartic) anharmonicity. The efficiency 
of these channels determines the width of the acousti- 
cally induced band gaps and the reflectivity spikes and 
increases with increasing acoustic intensity I ac ^- 

To model numerically the RAO effect in the LiNbOs 
filter wc calculate reflectivity from the acoustically modu- 
lated semiconductor as well as extinction of the THz light 
inside a LiNbOa crystal. The acoustic wave parametri- 
cally modulates the semiconductor through the macro- 
scopic polarization P(r, t) of TO phonons. This modu- 
lation can be expressed mathematically using the Hop- 
field models of polaritons, generalized to the case of a 
parametric perturbation by an external coherent acoustic 
wave. The material equation for the TO-phonon polar- 
ization is then extended to the form— 7 * 
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which includes both cubic and quartic terms of the cou- 
pling to the ultrasound wave, and in combination with 
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Maxwell's equation and boundary conditions allows us 
to calculate the macroscopic electric field E(r, £). The 
standard material parameters wto, 7tOi an d ur are the 
TO-phonon frequency, room temperature phenomenolog- 
ical TO-phonon damping, and the Rabi frequency of the 
Hopficld polariton, respectively. The background dielec- 
tric constant is given by E\>. There are two controlling pa- 
rameters for the RAO effect, the frequency of the acous- 
tic wave f2 ac = v ac /{2ir) = v s K, where K = |K| is the 
acoustic wave vector, and the strength of the coupling be- 
tween the acoustic wave and the TO-phonon field. This 
coupling strength is controlled by the acoustic intensity 

1/2 

I a c(r, t) and is given by m^I a c for cubic and m^I ac for 
quartic anharmonicity, the amplitudes of intrinsic cou- 
pling constants being 7723 = 1.52 x 10 7 s 1 / 2 g -1 / 2 and 
7714 = 1.64 s 2 g -1 in LiNb03. The coupling constants 
are derived from the potential of the anharmonic com- 
ponents added to the Hopfield Hamiltoniam 17 i 2( h 2£ l n 
this paper, the coupling to a plane bulk acoustic wave 
with a moderate intensity 7 ac (r, i) = I ac = 25kW/cm 2 is 
considered. The values of the phase shift <fi = ^3 — ^4 
between the complex anharmonic coupling parameters in 

Eq.Q]are not available in the literature. Despite the large 

1/2 

differences in coupling strength {hm^Iai = 5.02 meV 
while hm^I ac = 0.27 meV) and the fact that the quar- 
tic modulation alone leads to negligibly small changes 
in the LiNb03 reflection, the cubic and quartic modu- 
lations interfere with each other resulting in noticeable 
^-dependent effects which are also discussed below. 

The modulation of the medium due to either cubic or 
quartic anharmonicities has been studied by us in Refs.[l7l 
and[l||. In this Letter, we take into account both cubic 
and quartic modulations on an equal footing, including 
their interference. For collincar propagation of a THz 
light and an ultrasound wave, we calculate the Bragg 
scattering of the light by the acoustic wave, which is res- 
onantly enhanced by the TO-phonon component of the 
polariton in LiNb03. To do so, we use the numerical 
method developed in Refs.ll~7l. In particular, we solve nu- 
merically Eq. Q]in a semi-infinite LiNb03 region, using a 
superposition of partial polaritonic plane waves with all 
important Bragg harmonics taken into account. Because 
of the large strength of the mediated coupling we take 
into up to 400 harmonics. Applying Maxwcllian bound- 
ary conditions we calculate the amplitudes of the partial 
waves of the transmitted and reflected electric field. Full 
details of the calculation method can be found in Ref.ll8l. 

The total Bragg reflectivity which includes all impor- 
tant Bragg replicas, calculated for the extreme cases 
of constructive (4> = 0) and destructive interference 
(d> = 7r/2) is compared in Fig. [2] with the TO-phonon 
polariton reflectivity at I ac = (red dashed lines). The 
arrows indicate the peaks in the reflectivity caused by 
the gaps in the quasi-energy spectrum due to acousti- 
cally induced transitions between different TO-phonon 
polaritons. The peaks occur at the resonant Bragg con- 
dition lu = ujp i(nK/2) [wpoi(fc) is the TO-phonon polari- 
ton quasi-energy dispersion] , both in the lower polariton 
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FIG. 2. (Color online) (a) to (e) show the development the 
THz reflection of LiNbOg with varying acoustic frequency and 
different phase shift between the cubic and quartic anhar- 
monic parameters, = (black solid lines) and (f> — n/2 
(green dated lines). The red dashed line is the reflectiv- 
ity spectrum at 7 ac = 0. In both cases cjto = 31 meV, 
ur = 18.0 meV, 7x0 = 1.26 meV, and e b = 



(LP) and in the upper polariton (UP) branches, at the 
frequencies where co- and contra-propagating Bragg scat- 
tered polariton states are brought into resonancej 17 i 18 

The dramatic changes of the reflectivity peaks in a 
wide spectral range and, as suggested by the Bragg 
condition, their blue shift with increasing v ac demon- 
strated in Fig. [21 is the phenomenon which we propose 
to use for a widely tunable THz filter. The acoustic fre- 
quency provides an easy to manipulate parameter which 
allows for the tuning of the reflectivity. It can be clearly 
seen in Fig. [2] (a-d) that the light frequency range of 
1 < uj/(2ir) < 4 THz can be covered by the acoustic 
modulation of the LP branch with up to v ac ~ 1 GHz. 

The n-th peak in the reflection physically corresponds 
to stimulated emission of \n\ dressed acoustic phonons 
and to down-conversion of the incoming THz light. In 
the LP branch, such dressing with higher-order processes 
of re-emission and re-absorption of acoustic phonons by 
the TO-phonon polariton becomes more and more effi- 
cient with increasing acoustic frequency v ac , due to the 
flattening of the polariton dispersion near the bottom 
of the Rcstrahlen band and increasing of the TO-phonon 
fraction in the polariton states, which enhances their cou- 
pling to the acoustic wave. This results in widening of the 
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FIG. 3. (Color online) (a) The total Bragg reflectivity with 
the emphasis on the n — — 2 Bragg signal in LP branch cal- 
culated for i^ac = 250 MHz and different phase: = 0, tt/4, 
and 7r/2. (b) The cumulative cups-like peak in the reflectivity 
spectrum calculated for v ac = 2 GHz. 
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FIG. 4. (Color online) Extinction of the electric field at 
different acoustic frequencies, for = (black solid lines) 
and <j> — tv/2 (green dashed lines). 



lower peak with ^ ac and even more pronounced spectral 
changes in the Restrahlen band itself, see Fig. [2] (a-e). 

From Eq. Q] it is clear that the phase shift <\> will only 
affect the peak which is caused by both cubic and quartic 
anharmonicity. The cubic anharmonicity contributes to 
all Bragg signals (n = —1, n = —2, and so on), the quar- 
tic anharmonicity contributes to even n only. This is why 
the strongest interference is observed for n = —2 Bragg 
signal which is detailed in Fig. [3] (a) for different values of 
4>. The peak position remains practically unaffected by <p 
as it is determined within this acoustic frequency range 
solely by the resonant Bragg condition. For high acoustic 



frequency, shown in Fig. [3](b), the line width and height 
of the LP peak are not greatly effected by the phase. 
However, the peak position is strongly (/(-dependent be- 
cause it is no longer determined by the Bragg condition. 
Instead, the individual acoustically induced band gaps 
accumulate at the bottom edge of the Restrahlen band 
to produce a cusp-like feature seen in Fig. [3] (b) , which 
is discussed in detail in Ref.[l8l 

To quantify the transmission through the LiNbC>3 slab 
we plot in Fig. 4 the extinction which is calculated as an 
inverse decay length of the electric field inside the slab. 
The spectra show a high contrast in THz light filtering 
as well as its wide-range tunability that demonstrates a 
large figure of merit of the proposed device. 

In conclusion, our numerical modelling of the RAO ef- 
fect in LiNbC>3 shows that there is potential for a THz 
filter that can be continuously tuned between 0.5 THz 
and ~ 4 THz. This could be achieved through the ma- 
nipulation of the spectral position and height of peaks 
in reflectivity and extinction via the acousto-optical res- 
onant Bragg condition. 

The authors thank S. G. Tikhodccv and W. Langbcin 
for valuable discussions. The work was supported by the 
Royal Society (Grant JP0766306), EPSRC, and WIMCS. 

1 C. A. Schmuttenmaer, Chem. Rev. 104, 1759 (2004). 
2 J. C. Feurer, J. C. Vaughan, and K. A. Nelson, Science 299, 374 
(2003). 

3 S. Kojima, N. Tsumura, W. M. Takeda, and S. Nishizawa, Phys. 

Rev. B 67, 035102 (2003). 
4 P. Y. Han and X.-C. Zhang, Meas. Sci. amd Technol. 12, 1747 

(2001). 

5 N. S. Stoyanov, D. W. Ward, T. Feurer, and K. A. Nelson, 

Nature Mat. 1, 95 (2002). 
6 P. H. Siegel, IEEE T. Microw. Theory 50, 910 (2002). 
7 B. Fischer, M. Hoffmann, H. Helm, G. Modjesch, and P. U. 

Jepsen, Semicond. Sci. Technol. 20, S246 (2005). 
8 J. F. Federici, B. Brian Schulkin, F. Huang, D. Gary, R. Barat, 

F. Oliveira, and D. Zimdars, Semicond. Sci. Technol. 20, S266 

(2005). 

9 I. H. Libon, S. Baumgartncr, M. Hemplc, N. E. Hecker, J. Feld- 

mann, M. Kock, and P. Dawson, Opt. Lett. 33, 1401 (2008). 
10 C.-L. Pan, C.-F. Hsieh, R.-P. Pan, M. Tanaka, F. Miyamaru, 

M. Tani, and M. Hangyo, Opt. Express. 13, 113921 (2005). 
n C.-Y. Chen, C.-L. Pan, C.-F. Hsieh, Y.-F. Lin, and R.-P. Pan, 

Appl. Phys. Lett. 88, 101107 (2006). 
12 I.-C. Ho, C.-L. Pan, C.-F. Hsieh, and R.-P. Pan, Opt. Lett. 33, 

1401 (2008). 

13 T. D. Drysdale, I. S. Gregory, C. Baker, E. H. Linfield, W. R. 

Tribe, and D. R. S. Cumming, Appl. Phys. Lett. 85, 5173 (2004). 
14 A. L. Ivanov and P. Littlewood, Phys. Rev. Lett. 87, 136403 

(2001). 

15 A. L. Ivanov, Phys. Status. Solidi A 202, 2657 (2005). 

16 A. L. Ivanov, in Problems of Condensed Matter Physics, edited 

by A. L. Ivanov and S. G. Tikhodeev (Oxford University Press, 

Oxford, 2008) pp. 301-322. 
17 R. H. Poolman, E. A. Muljarov, and A. L. Ivanov, Phys. Rev. 

B 81, 245208 (2010). 
18 E. A. Muljarov, R. H. Poolman, and A. L. Ivanov, (2011), in 

print. 

19 J. J. Hopfield and D. G. Thomas, Phys. Rev. 132, 563 (1963). 
20 V. Romero-Rochm, R. M. Kochl, C. J. Brennan, and K. A. 

Nelson, J. Chem. Phys. Ill, 3559 (1999). 
21 S. C. Abrahams, in Properties of Lithium Niobate, edited by 

K. K. Wong (The Institute of Electrical Engineers, 2002). 
22 J. Menendez and M. Cardona, Phys. Rev. B 29, 2051 (1984). 



